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Magnetic materials and nanostructures based on carbon offer unique opportunities for future 
technological applications such as spintronics. This article reviews graphene-derived systems in 
which magnetic correlations emerge as a result of reduced dimensions, disorder and other possi- 
ble scenarios. In particular, zero-dimensional graphene nanofragments, one-dimensional graphene 
nanoribbons, and defect-induced magnetism in graphene and graphite are covered. Possible phys- 
ical mechanisms of the emergence of magnetism in these systems are illustrated with the help of 
computational examples based on simple model Hamiltonians. In addition, this review covers spin 
transport properties, proposed designs of graphenc-based spintronic devices, magnetic ordering at 
finite temperatures as well as the most recent experimental achievements. 
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I. INTRODUCTION 

Magnetic materials are essential for modern technol- 
ogy. All presently used magnetic materials involve the 
elements belonging to either the d- or the /-block of the 
periodic table. For instance, among the periodic table el- 
ements only the late transition metals Fe, Co and Ni, are 
ferromagnets at room temperature. Magnetic ordering in 
these transition metals originate from the partially filled 
d-electron bands. However, magnetism is not common 
for the light p-block elements belonging to the second 
period of the periodic table, even despite the fact that 



carbon is able to form very diverse and complex molec- 
ular structures. In principle, such materials may pos- 
sess a number of attractive properties such as low den- 
sity, biocompatibility, plasticity, and many others, which 
stimulates the search for light -element based magnetism 
(jMakarova and Palaciol 120061 ). 

The field of light-element magnetism and, in particu- 
lar, of carbon-based magnetism is currently gaining in- 
creasing importance because of the following two reasons. 
Firstly, the field of carbon-based magnetism has always 
been a very controversial area of research which suffered 
from the poor reproducibility of experimental results. 
However, the situation seems to be improved over the 
last few years. Several examples of magnetism in carbon- 
based materials continue to be reliably reproduced by 
different research groups. The second driving force is the 
first isolation of graphene, a truly two-dimensional form 
of carbon which has a ttracted enormous atte ntion in sci- 
ence and technology dNovoselov et al I I2004T) . Graphene 
has a fairly simple honeycomb atomic structure, but 
rather unique electronic structure with linear band dis- 
persion at the Fermi level (see Fig. rjj which is largely 
responsible for many novel phy sical phenomena observed 
in this material (for review see (ICastro Neto et all 120091: 
iGeim and Novoselovi l2007t iKatsnelsonl . 120071) ). Impor- 
tantly. graphene can also be considered as a unifying 
concept for understanding a broad class of sp 2 carbon 
materials which includes polycyclic aromatic molecules, 
fullerenes, carbon nanotubes, and graphite as well as 
their further modifications obtained by patterning, chem- 
ical treatment, implantation of defects, impurities, etc. 

While ideal graphene is nonmagnetic itself, many its 
derivative materials and nanostructures, both realized in 
practice and considered theoretically, show various sce- 
narios of magnetism. The magnetic graphene nanostruc- 
tures are particularly promising for applications in the 
field of spintronics, a very probable future step in the evo- 
lution of electronics industry, which promises information 
storage, processing and commu nication at faster speeds 
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and lo wer energy cons u mptio n (lAwschalom and Flatte 



120071 : IChappert et al. |, 120071 : iFertl. l2008t IWolf et al 
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FIG. 1 (a) Two-dimensional crystalline lattice of graphene. 
The shaded area denotes the unit cell of graphene contain- 
ing two carbon atoms which belong to the two sublattices of 
graphene, A (empty circle) and B (filled circle). The two 
high-symmetry directions of graphene lattice, armchair and 
zigzag, are highlighted, (b) Band structure of graphene ob- 
tained by means of first-principles calculations. The bands 
are labeled according to their symmetry. The 7r-symmetry 
bands responsible for the low-energy electronic properties of 
graphene are highlighted. The zero energy corresponds to the 
Fermi level, (c) The low-energy part of the band structure 
of graphene involves two inequivalent 'Dirac cone' features 
in the corners (points K and K') of the hexagonal Brillouin 
zone. 



and decoherence of electron sp i ns (iFischer et all 120091 : 
iTrauzettel et all 120071 : lYazyevl l2008aj). This property 
makes carbon nanomaterials promising for transport of 
spin-polarized currents and for spin-based quantum in- 
formation processing. 

This review provides a brief introduction into the cur- 
rent state of the field of magnetic materials and nanos- 
tructures based on graphene. The possible scenarios for 
the onset of magnetism in graphene nanostructures are 
illustrated by means of a simple theoretical model based 
on the mean-field Hubbard Hamiltonian. The next sec- 
tion briefly reviews the landmark experimental reports 
in the field. Then follows the introduction of the the- 
oretical model and its specific consequences for describ- 
ing the electronic structure and magnetic properties of 
graphene-based materials. The main part of this arti- 
cle applies the theoretical model described and reviews 
both theoretical and experimental published works on 
magnetic graphene systems classified according to their 
dimensionality: (I) finite graphene nanofragmcnts, (2) 
one-dimensional graphene edges and graphene nanorib- 
bons, and (3) two-dimensional graphene and graphite 
with magnetism induced by the presence of point de- 
fects. Of these three classes, magnetic graphene edges 
and nanoribbons will be covered in more detail since 
these one-dimensional objects continue to receive special 
attention in the scientific community. Future perspec- 
tives of the field arc outlined in the last section. 



II. BRIEF OVERVIEW OF EXPERIMENTAL PROGRESS 

The first reproducible experimental reports of 
magnetism in p-block compounds were published 
in 1991 when magnetic ordering was observed 
in crystalli ne p-nitropheny l nitronyl nitroxide (p- 



NPNN) (jTakahashi et al 
Il99lh and 



Il99ll iTamura et al\ . 
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a charge transfer complex of C 
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l200ll ). While traditional electronics exploits only the 
charge of electron, spintronics will also make use of its 
spin degree of freedom. For the field of spintronics 
graphene can offer a possibility of tuning its spin trans- 
port properties by means of various applied stimuli. For 
instance, it was suggested that half-metallicity of zigzag 
graphene nanoribbons can b e triggered by external elec- 
tric fields ( Son et al l l2006bl) . If realized in practice, this 
would allow for efficient electric control of spin trans- 
port, a very desirable effect which is hard to achieve 
using other materials. In addition, materials based 
sp-elements arc expect ed to have high magnitudes o f 
the spin- wave stiffness ( Edwards and Katsnelsonl . 120061 ) 



and, thus, nanostructures made of these elements would 
possess higher Curie temperatures o r spin correlation 
lengths ( Yazvev and Katsnelsonl |2008| ) . Materials based 
on light elements also display weak spin-orbit and hyper- 
fine couplings which are the main channels of relaxation 
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FIG. 2 Chemical structures of (a) p-NPNN and (b) the 
TDAE-Ceo charge transfer complex. Molecular crystals of 
these organic compounds exhibit ferromagnetic ordering with 
Curie temperatures of 0.6 K and 16 K, respectively. Each 
molecular radical unit bears a single uncompensated electron 
spin. In chemical notation the presence of aunpaired electron 
is usually denoted by a thick dot. 
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and tetrakis ( dimeth ylamino)cthylene (TDAE) 

(jAllemand et all Il99ll) . Molecular structures of 



these organic materials made of light elements only (C, 
H, N and O) are shown in Figured In their molecular 
crystals the uncompensated electron spins are localized 
on weakly coupled molecular units. Because of the 
weak coupling between electron spins, the long-range 
magnetic order is realized only at low temperatures. The 
two organic materials mentioned above, p-NPNN and 
TDAE-Cgo, are characterized by Curie temperatures of 
0.6 K and 16 K, respectively. Since 1991 a large number 
of other organic magnetic materials have been examined. 
In all cases the temperatures below which long-range 
magnetic order is established (Curie temperatures, Tq, 
and Neel temperatures, Tn, in the case of ferromagnetic 
and antiferromagnetic orderings, respectively) were 
much lower than room temperature, which renders such 
materials useless for practical applications. 

The next milestone experiment was reported ten 
years later when ferromagnetism with Tq w 500 K 
was observed in rhomboh edral Ceo under high pres- 
sure ( Makarova et al 1. l200lh . This observation, however, 
demonstrates very well the controversial character of the 
field. Five years later several authors retracted the orig- 
inal publication since the measured content of magnetic 
impurities was shown to be close to the amount needed 
to explain the observe d magnetization of the samples 
(|Makarova et all l2006lh [n addition, the measured Tc 
was found to be very similar to the one of cementite 
FC3C. The question of possible high-temperature mag- 
netic ordering in Ceo-based materials remains open. 

Two years later room-temperature ferromagnetism 
was observed in highly oriented pyrolytic graphite 
(HO PG) irradiated with hi gh-energy (2.25 MeV) pro- 
tons (jEsauinazi et all 120031 ). Figure [3] shows the mag- 
netization loop for a proton-irradiated sample compared 
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FIG. 3 Magnetic moment of proton irradiated (filled circles) 
and untreated (empty squares) graphite as a function applied 
magnetic field as measured by Esquinazi et al. at T — 300 K. 
The weak hyster e sis loo p can be recognized. Reprinted from 
(|Esquinazi et all I2003T ). Copyright 2003 by the American 
Physical Society. 



to untreated HOPG. Further experimental investigations 
revealed that the magnetic order in proton-bombarded 
graphi te has two-dimensional , that i s, graphene-like char- 
acter ( Barzola-Quiquia et all I2007I ) and originates from 
the carbon 7r-clcctron s ystem rather than f rom possi- 
ble d-element impurities ( Ohldag et all l2~007t ) . Interest- 
ingly, it was shown that the chemical nature of the high- 
energy particles plays a crucial role in producing mag- 
netic ordering. While proton irradiation leads to the on- 
set of ferromagnetism i n irradiat ed samples, both helium 
(lEsquinazi et al\ , I2003I ) and iron ( Barzola-Quiquia et all 
I2008D ions show no clear effect. On the other hand, the 
implantation of carbon ion s was also repor ted to induce 
ferromagnetism in HOPG (|Xia et aUl2008D . 

A number of reports have also pointed out 
that even untreated graphite exhibit f e rroma gnetism 
(|Esquinazi et all I2OO2I ; iKopelevich et all l2000lh Very 



recently, by using a combination of scanning probe tech- 
niques and magnetization measurements, Cervenka and 
co-authors have shown that the intrinsic ferromagnetism 
of graphite is related to the presence of grain bound- 
aries which c an be considered as 2D periodic networks of 
point defects ( Cervenka et «Ll . l2009l ). Room-temperature 
magnetic hysteresis has also been reported for graphene 
samples produced in bul k quantities from graphite using 
the ch emical approaches ( Matte et aLl . l200"9tlWang et all 
2009blh 



III. BASIC COMPUTATIONAL APPROACHES 

A. Model Hamiltonians 

The vast majority of computational studies of mag- 
netic carbon nanostructures are currently performed us- 
ing first-principles electronic structure methods based 
on density functional theory (DFT). This approach is 
now implemented in a large number of public computer 
codes and well described i n a variety of graduate-level 
textbooks (for instance, see dKoch and Holthausenl . l2002t 
iMartinl . |2004 IMarx and Hutterl . l2009l) h For pedagogical 
purposes, a simpler approach based on model Hamilto- 
nians will be adopted in this review article. Moreover, 
it will be demonstrated below that such simplified mod- 
els very often allow deeper understanding of the results 
obtained. 

A simple model which is widely used for studying mag- 
netic effects in sp 2 carbon materials is the one-orbital 



mean-field Hubbard model. This model considers only 
the 7r-symmetry electronic states which are formed by 
the unhybridized p z atomic orbitals of sp 2 carbon atoms. 
As shown in Figure [ljb) the low-energy electronic states 
have 7r-symmetry and thus play the dominant role in 
the properties of graphene systems. The Hubbard model 
Hamiltonian can be partitioned into two parts, 

H = H +H'. (1) 

The first term is the nearest-neighbor tight-binding 
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in which the operators Ci a and c icr annihilate and create 
an electron with spin a at site i, respectively. The nota- 
tion (-, •) stands for the pairs of nearest-neighbor atoms; 
'h.c' is the Hcrmitian conjugate counterpart. The well 
established hopping integral t w 2.7 eV defines the en- 
ergy scale of the Hamiltonian. This physical model is 
equivalent to the Hiickel method familiar to chemists. 
From the computational point of view, the Hamiltonian 
matrix is determined solely by the atomic structure: the 
off-diagonal matrix elements and are set to 

—t when atoms i and j are covalently bonded, and to 
otherwise. In a neutral graphene system each sp 2 car- 
bon atom contributes one p z orbital and one 7r electron. 
The 7r-electron system is thus called half-filled. The spec- 
trum of the eigenvalues of tight-binding Hamiltonian ma- 
trix exhibit electron- hole symmetry, i.e. it is symmetric 
with respect to zero energy. In other words, in a neutral 
graphene system for each eigenvalue e < correspond- 
ing to an occupied (bonding) state, there is an unoccu- 
pied (anti-bonding) state with e* = — e. The states with 
e = are called zero-energy states (also referred to as 
non-bonding or midgap states). 

The nearest-neighbor tight-binding model has proved 
to describe accurately the electronic structure of 
graphene, carbon nanotubes and other non-magnetic sp 2 
carbon materials. However, electron-electron interac- 
tions have to be introduced in some form in order to 
describe the onset of magnetism. Within the Hubbard 
model these interactions are introduced through the re- 
pulsive on-site Coulomb interaction 



(3) 



where = ct c, CT is the spin- resolved electron density 
at site i; the parameter U > defines the magnitude of 
the on-site Coulomb repulsion. This model considers only 
the short-range Coulomb repulsion, that is, two electrons 
interact only if they occupy the p z atomic orbital of the 
same atom. Despite its apparent simplicity, this term is 
no longer trivial from the computational point of view. 
The mean-field approximation 

^mf = U {n^{n iir ) + (nit)nu - (n^}{n ti )) , (4) 



allows us to overcome this difficulty. Here, a spin-up 
electron at site i interacts with the average spin-down 
electron population (n^) at the same site, and vice 
versa. This mean-field model represe nts a variation of the 
unres tricted Hartree-Fock method ( Szabo and Ostlundl . 
Il982h . From the computational point of view, the 
electron-electron interaction term affects only the diag- 
onal elements of the Hamiltonian matrix. The diago- 
nal elements of the spin-up and spin-down blocks now 



The problem can be solved self-consistently starting from 
some initial values of (rii a ) which can be chosen ran- 
domly. However, one has to keep in mind that in certain 
cases the broken-symmetry (antiferromagnetic) solutions 
can be obtained o nly if the initia l gues s breaks the spin- 
spatial symmetry ( Yazvev et al\ . 12008ft . The process of 
(I) calculation of the matrix elements of the Hamiltonian 
matrix, (2) its diagonalization and (3) the computation 
of updated spin densities is then repeated iteratively until 
all values of (rii a ) are converged. The final self-consistent 
solution provides the spin densities 



Mi 



{mi) 



(5) 



at each atom i and the total spin of the system S = 
Mj. For a given graphene structure both local and 
total spins (magnetic moments) depend exclusively on 
the dimensionless parameter U /t. 

After the model has been introduced, the following 
three critical questions can be asked. (I) Is the one- 
orbital approximation accurate enough compared to the 
methods considering all electrons? (2) Which value of 
the empirical parameter U /t should be used? (3) Is the 
mean-field approximation justified for graphene based 
systems? 

(1) It has been shown that the results of mean- 
field Hubbard model calculations correspond closely 
to the ones obtained using first-principles meth- 
ods if the parameter U/t is chosen appropriately 
(iFernandez-Rossier and Palaciost l2007t iGunlvck e et all 
120071 ; iPisani et all 120071 1 The first-principles meth- 
ods either treat all electrons on equal footing or disre- 
gard the localized atomic core states which are not im- 
portant in most cases. One notable exception is the 
calculation of hyperfine interactions. In this case the 
spin polarization of the Is atomic core states of car- 
bon atoms has significant contribution to the Fermi con- 
tact hyperfine couplings ( Yazvevl . l2008at lYazvev et all 
12005ft . Otherwise, the results of DFT calculations per- 
formed using a generalizcd-gradient-approximation fam- 
ily exchange-correlation functional are best reproduced 
when U/t w 1.3. The results of the local-spin-density ap- 
proximation c alcula tions are best fitted using U/t ps 0.9 
(jPisani et aTl . l2007ft . 

(2) Ideally, the empirical parameter U/t must be es- 
timated using experimental knowledge. Unfortunately, 
there are no direct experiments performed on mag- 
netic graphene systems which would allow to estimate 
U/t. Magnetic resonance studies of neutral soliton 
states in irans-polyacetylene, a one-dimensional sp 2 car- 
bon polymer which can be viewed as a minimum- 
width zigzag graphene n anoribbon, give the range of 
values U ~ 3.0—3.5 e V ( Kuroda and Shirakawal Il987t 
iThomann et al. I, 11985ft . This interval corresponds to 
U/t ~ 1.1—1.3 which also makes us confident in the 
results of the gencralized-gradicnt-approximation DFT 
calculations. Increasing U /t leads to the enhancement 
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of magnetic moments. The range of meaningful mag- 
nitudes is limited by U/t w 2.23 above which the 
ideal graphene undergoes a Mott-Hubbard transition 
i nto an antiferromagnetic ally ordered insulating state 
(jSorella and Tosattil . 119921 ). In the computational exam- 
ples considered below a value of U/t = 1.2 will be used. 

(3) This question is the most difficult to answer. A 
comparison of the mean-field results with the ones ob- 
tained using exact diagonalization and quantum Monte 
Carlo simulations illustrates the validity of this approx- 
imati on for the relevant values of U/t (jFeldner et all 
2009). Magnetic graphene materials and nanostructures 
need not be considered as strongly correlated systems. 



instability associated with the presence of low-energy 
electrons in the system. Other mechanisms, such as the 
Peierls distortion, were show n to be inefficient in the case 
of graphene nanostructures ( Pisani et all l2007t ) . 

Although the benzenoid graph theory is able to predict 
the occurrence of zero-energy states, it is not clear how 
the electron spins align in these states. The complemen- 
tary knowledge is supplied by Lieb's theorem (jLiebl . [l989h 
which determines the total spin of a bipartite system 
described by the Hubbard model. This theorem states 
that in the case of repulsive electron-electron interactions 
(U > 0), a bipartite system at half- filling has the ground 
state characterized by the total spin 



B. Counting rules 

There are two important consequences coming from 
the model Hamiltonians we have introduced. The hon- 
eycomb lattice of graphene is a bipartite lattice. That is, 
it can be partitioned into two mutually interconnected 
sublatticcs A and B (see Fig. QJa) ) . Each atom belong- 
ing to sublattice A is connected to the atoms in sublatticc 
B only, and vice versa. Moreover, the graphene systems 
whose faces are hexagons are called benzenoid (or hon- 
eycomb) systems. Carbon atoms in such systems have 
either three or two nearest neighbors. The class of ben- 
zenoid systems is a subclass of bipartite systems. 

The spectrum of the tight-binding Hamiltonian of a 
honeycomb system can be analyzed using a mathemat- 
icall y rigorous approach of the benzenoid graph the- 
ory ( Faitlowicz et all , |2005f). An important result for 
us is that this theory is able to predict the number of 
zero-energy states of the nearest-neighbor tight-binding 
Hamiltonian in a 'counting rule' fashion. The number of 
such states is equal to the graph's nullity 



S=-\N A -N B \, 



(8) 



77 = 2a - N, 



(6) 



where N is the total number of sites and a is the maxi- 
mum possible number of non-adjacent sites, i.e. the sites 
which are not the nearest neighbors to each other. The 
onset of magnetism in the system is determined by the 
so-called Stoner criterion which refers to the competi- 
tion of the exchange energy gain and the kinetic-energy 
pena lty associate d with the spin-polarization of the sys- 
tem (lMohnl . l2003h . The gain in exchange energy is due to 
the exchange splitti ng of the electronic st ates subjected 
to spin-polarization ( Palacios et all 120081 ) 



At 



U 



(7) 



where ^ nf is the inverse participation ratio, a mea- 
sure of the degree of localization of the corresponding 
electronic state. The kinetic-energy penalty is propor- 
tional to the energy of this state. Thus, the zero-energy 
states undergo spin-polarization at any U > irrespec- 
tive of their degree of localization. One can view spin- 
polarization as one of the mechanisms for escaping an 



where Na and Nb are the numbers of sites in sublat- 
tices A and B, respectively. The ground state is unique 
and the theorem holds in all dimensions without the ne- 
cessity of a periodic lattice structure. Importantly, the 
two counting rules are linked by the following relation, 
V> \N A -N B \. 

In the following section, the application of these two 
simple counting rules will be demonstrated on small 
graphene fragments and compared to the results of nu- 
merical calculations. 



IV. FINITE GRAPHENE FRAGMENTS - A SIMPLE 
ILLUSTRATION 

Let us first try to understand the origin of mag- 
netism in finite graphene fragments (also referred to as 
nanoflakes, nanoislands or nanodisks) as a function of 
their shape and size. Three simple examples of nanome- 
ter sized graphene fragments are shown in Figured From 
the point of view of single-orbital physical models, only 
the connectivity of the 7r-electron conjugation network 
is important. Such 7r-systems may constitute only small 
parts of more complex molecules or bulk materials. In 
simplest case, the 7r-electron networks shown in Figure 0] 
can be realized in the corresponding all-benzenoid poly- 
cyclic aromatic hydrocarbon (PAH) molecules with the 
edges of the fragments being passivatcd by hydrogen 
atoms. Each carbon atom at the edge of the fragment 
is bonded to one hydrogen atom such that all carbon 
atoms are sp 2 -hybridized. Current progress in synthe- 
sizing such molecules and u nderstanding the ir properties 
has recently been reviewed ( Wu et aZl . l2007|) . 

The hexagonal graphene fragment shown in Fig- 
ure SJa) is thus equivalent to the coronene molecule. For 
this fragment, the number of sites belonging to the two 
sublattices is equal, Na = Nb = 12. The number of non- 
adjacent sites is maximized when all atoms belonging to 
either of the two sublattices are selected, i.e. a = 12. 
Thus, both the number of zero-energy states rj and the 
total spin S are zero. The tight-binding model predicts 
a wide band gap of 1.08i ~ 3.0 eV for this graphene 
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molecule. As expected, the mean-field Hubbard model 
solution for this fragment does not reveal any magnetism. 







0.5 
band filling 

FIG. 4 Atomic structures and tight-binding energy spectra 
of three graphene fragments: (a) coronene, (b) triangulane, 
and (c) a bowtie-shaped fragment ("Clar's goblet"). Non- 
adjacent sites are labeled by circles. Empty and filled circles 
correspond to sublattice A and sublattice B, correspondingly. 
Tight-binding energies are plotted as a function of band fill- 
ing. Dashed line corresponds to the energy spectrum of ideal 
graphene. 

The second graphene fragment shown in Figure Eth) 
has triangular shape. It is not surprising that the cor- 
responding hypothetical PAH molecule is called triangu- 
lane. Unlike coronene, the two sublatticcs of this trian- 
gular fragment are no longer equivalent: Na = 12 and 
Nb = 10. The unique choice maximizing the number of 
non-adjacent sites is achieved by selecting the atoms be- 
longing to the dominant sublattice A, i.e. a = Na = 12. 
Thus, the benzenoid graph theory predicts the presence 
of two zero-energy states on sublattice A. Lieb's theorem 
predicts the S = 1 (spin-triplet) ground state or, cquiv- 
alcntly, a magnetic moment of 2 /ie per molecule. The 
two low-energy electrons populate a pair of zero-energy 



states according to Hund's rule, that is, their spins are 
oriented parallel to each other. The mean-field Hubbard 
model results for this system at half-filling are shown 
in Figure [5ja) . One can see that spin-polarization lifts 
the degeneracy of the zero-energy electronic states and 
opens an energy gap As = 0.30t sa 0.8 eV. The system 
is stabilized by spin-polarization. Most of the spin-up 
electron density localized on the atoms in sublattice A 
(see Fig.^a)) originates from the two electrons populat- 
ing the non-bonding states. However, one can notice an 
appreciable amount of spin-down density on the atoms 
in sublattice B which is compensated by an equivalent 
contribution of the spin-up density in sublattice A. The 
occurrence of the induced magnetic moments is a man- 
ifestation of the spin-polarization effect which is related 
to the exchange interaction of the fully populated states 
with the two unpaired electrons. 

The third bowtie-shaped graphene molecule shown in 
Figure Etc) is composed of two triangulane fragments 
sharing one hexagon. For this system Lieb's theorem 
predicts the spin-singlet ground state (Na = Nb = 19). 
However, the choice of the set of atoms which maxi- 




0.75 



0.75 



FIG. 5 Local magnetic moments and spin-resolved energy lev- 
els obtained through the mean-field Hubbard model calcula- 
tions for (a) triangulane and (b) the bowtie-shaped graphene 
fragment (U/t — 1.2). Area of each circle is proportional to 
the magnitude of the local magnetic moment at each atom. 
Filled (red) and empty (blue) circles correspond to spin-up 
and spin-down densities. Energy levels energies are plotted 
as a function of band filling. Dashed line corresponds to the 
energy spectrum of ideal graphene. Red and blue levels corre- 
spond spin-up and spin-down channels, respectively. In the 
case of bowtie fragment the energies in the two channels are 
identical (shown in magenta). 
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mizes the number non-adjacent sites is less evident in 
this case. Figure SJc) shows such a selection (a = 20) 
which involves the atoms belonging to both sublattice A 
and sublattice B in the left and right parts of the struc- 
ture. These atoms are marked differently in the figure. 
Hence, there are 77 = 2x 20 — 38 = 2 zero-energy states 
as confirmed by the tight-binding calculation. The zero- 
energy states are spatiall y segregated in the t wo triangu- 
lar parts of the molecule dWang et all l2009al) . To satisfy 
the spin-singlet ground state, the two zero-energy states 
have to be populated by two electrons with oppositely 
oriented spins. In other words, the ground electronic 
configuration breaks spin-spatial symmetry and exhibits 
antiferromagnetic ordering. This result can be verified 
by mean-field Hubbard model calculations as shown in 
Figure E^b). It can be argued that this example violates 
Hund's rule. However, one has to keep in mind that 
each of the two non-bonding states is localized within 
one of the graphene sublattices. That is, there are two 
electronic sub-bands, each populated by electrons ac- 
cording to Hund's rule. The coupling between the elec- 
tron spins in these two sub-bands is antiferromagnetic 
due to the sup erexchange mechanism ( Anderson] . Il950t 
lKramersl . fl93l . 

The two counting rules can be applied to larger 
graphene fragments. It was shown that the total spin 
of triangular fragments with edges cut along th e zigzag 
direction scales linearly with fragment size (lEzawa , 



20071; iFernandez-Rossier and Palaciosl 120071 1 Wang etal. , 



2008al ) as illustrated in Figure lol The average magnetic 
moment per carbon atom thus decays with increasing the 
system size. The evolution of magnetic properties with 
increasing size for hexagonal fragments with edges cut 
along the same zigzag direction is less trivial. It has been 
shown theoretically that above some critical size the sys- 
tem undergoes a transi tion into a broken-symmetry an- 
tiferr omagnetic state ( Fernandez- Rossier and Palaciosl 
I2007D . The critical size itself depends strongly on the 
value of U/t. However, it is easier to explain the ori- 
gin of this behavior in large systems from the standpoint 
of edge magnetism, which will be explained in the next 
section. 

Finally, a few words have to be said about the possi- 
bility of realizing in practice the magnetic graphene frag- 
ments we have discussed. It is expected that such mag- 
netic systems are more reactive than the non-magnetic 
polyaromatic molecules. Although triangulane itself has 






FIG. 7 Chemi c al de rivativ es of triangulane s ynthesized by 
(|Allinson et all Il995f ) and (|lnoue et all 1200 if ) . Their spin- 
triplet ground state has been verified by means of electron 
spin resonance measurements. 



never been isolated, successful synthesis of its chem- 
i cal derivatives shown in Figure IT ] has been reported 
(lAllinson et all Il995l ; llnoue et all\200l\ ). The spin- 



119951 : llnoue et azFl200lh . 

triplet ground state of these chemical compounds was 
verified by the electron spin resonance measurements. 
In principle, this example can be considered as an in- 
direct proof of edge magnetism in graphene systems, at 
least in finite fragments produced by means of the chem- 
ical bottom-up approach. The synthesized triangulane 
derivatives are reactive molecules, but nevertheless can 
be handled in common organic solvents and stored for 
many months at room temperature pro vided the solu- 
tion i s isolated from atmospheric oxygen ( Allinson et all 
Il995f) . Larger magnetic triangular molecules have not 
been synthesized so far. The PAH molecule correspond- 
ing to the considered bowtie fragment was hypoth esized 
by E ric Clar and named "Oar's goblet" after him ( Clarl . 
Il972f). Attempts to synt hesize this molecule have failed 
(IClar and Mackavl . Il972fh 

The examples shown above illustrate how three dif- 
ferent magnetic scenarios can be realized in very sim- 
ple finite graphene systems. These examples also pro- 
vide a way for designing nanostructurcs with predefined 
magnetic interactions, a highly useful tool for developing 
novel spintronic devices. The value of this approach has 
already been demonstrated by the proposal of reconfig- 
urable spintronic logic gates exploiting the strong anti- 
ferromagne tic couplings in the bowtie-shaped graphene 
fragments (|Wang et all l2009al) . Several devices for con- 
trolling spin-currents based on triagular graphene frag - 
ments have also been described recently ( Ezawai l2009bl) . 
A number of other intriguing properties predicted for 
zero-dimensional graphene fragments are potentially in- 
terest ing from bo t h fundamental and applie d points of 
view (lAkola et aZ. l. l2008l : lEzawaLl2008ll2009alh 
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V. ZIGZAG EDGES AND NANORIBBONS 
A. Physical mechanism of edge magnetism 



FIG. 6 Evolution of the total spin of triangular graphene 
fragments with size. 



As one moves on towards larger graphene fragments 
or infinite systems, the application of counting rules be- 
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unit cell 



FIG. 8 Atomic structures and tight-binding band structures 
of (a) armchair and (b) zigzag graphene nanoribbons. Unit 
cells of the graphene nanoribbons are shown. 



zigzag edge) as shown in Figure [S^b). Strictly speaking, 
the flat band does not correspond to zero-energy states, 
but rather to the states with energies the approach zero 
with increasing nanoribbon width. The low-energy states 
are localized at the edge and decay quickly in the bulk. 

High density of low-energy electronic states suggests a 
possibility of magnetic ordering. Indeed, the mean-field 
Hubbard model solution for this system reveals mag- 
netic moments localized at the edges as shown in Fig- 
ure [9ja). The localized magnetic moments display fer- 
romagnetic ordering along the zigzag edge while the mu- 
tual orientation of the magneti c moments localized at the 
opposite edges is antiparallel dFernandez-Rossierl 120081 : 
iFuiita et a/fll996l;ISon et aUl2006bl ). Thus, the net mag- 
netic moment of a zigzag nanoribbon is zero in agreement 
with Lieb's theorem (Na = Nb)- The band structure 
corresponding to the mean-field Hubbard model solution 
is compared to the tight-binding band structure in Fig- 
ure EJb). The introduced electron-electron interactions 
open a band gap across the whole flat-band segment 




comes impractical. An alternative approach considers 
the effects of edges of graphene nanostructures which 
can be conveniently modeled using one-dimensional pe- 
riodic strips of graphene. Such models are commonly 
referred to as graphene nanoribbons. There are two 
high-symmetry crystallographic directions in graphene, 
armchair and zigzag, as shown in Figure QJa) . Cutting 
graphene nanoribbons along these directions produces 
armchair and zigzag nanoribbons, respectively (Fig. [8|). 

The band structures of armchair and zigzag nanorib- 
bons are remarkably different. Figure Eta) shows the 
tight-binding band structure of a ^1.5 nm wide arm- 
chair nanoribbon. For this particular armchair graphene 
nanoribbon, introducing a pair of parallel armchair edges 
opens a gap of 0.26i. The nearest-neighbor tight- 



ing behavior for armchair nanoribbons (IBarone et at. 


2006; Brev et al. 


20071 Ezawal 2006 


;lNakada et al. 1996 


Peres et al. 2006 


: Son et all 12006a 


), and the two situa- 



tions alternate as the nanoribbon's width increases. The 
band gap of semiconducting graphene nanoribbons de- 
creases with increasing width. In the case of metallic 
nanoribbons two bands cross the Fermi level at the T 
point. No magnetic ordering is predicted in this case. 

Within the same model all zigzag graphene nanorib- 
bons are metallic and feature a flat band extending 
over one-third of the one-dimensional Brillouin zone at 
k G (2ir / 3a; ir / a) (a — 0.25 nm is the unit cell of the 




FIG. 9 (a) Local magnetic moments in a zigzag graphene 
nanoribbon calculated using the mean-field Hubbard model 
(U/t = 1.2). Area of each circle is proportional to the magni- 
tude of the local magnetic moment at each atom. Filled (red) 
and empty (blue) circles correspond to spin-up and spin-down 
densities, respectively, (b) Mean-field Hubbard-model band 
structure (solid magenta lines) compared to the tight-binding 
band structure (dashed lines) for the solution shown in panel 
(a). The band structures for spin-up and spin-down electrons 
are equivalent, (c) Mean-field Hubbard-model band structure 
for the same graphene nanoribbon with the ferromagnetic mu- 
tual orientation of the edge spins. The band structures for 
the majority-spin electrons and the minority-spin electrons 
are shown as red and blue lines, respectively. 
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turning the system into a semiconductor (Ag = 0.20t 
at U/t = 1.2). The spin-polarization almost does not 
affect the electronic states at higher energies. The band 
structures for the two spin channels are equivalent, but 
spin-spatial symmetry is broken. 

The coupling between the magnetic edges can be as- 
cribed to the superexchange mechanism as in the case 
of bowtie graphene fragment considered in the previous 
section. The magnitude of the antiferromagnetic cou- 
pling shows a w~ 2 depend e nce a s a function of nanorib- 
bon width w (|Jung et all 120091 ). The intercdgc mag- 
netic coupling strength of ~25 meV has been calcu- 
lated h^m_fixs_^principles for a ~1.5 nm wide nanorib- 
bon (IPisani et al ., 2007). Unlike the antiferromagnetic 
ground state, a zigzag graphene nanoribbon with ferro- 
magnetic interedge orientation is a metal with two bands 
crossing the Fermi level at k « 27r/3a (Fig. (Hie)). The 
possibility of switching between the two states was ex- 
ploited in a proposal of a gr aphene-based magnetic sensor 
( Munoz-Roias et all . l2009l ) described in the next section. 
The coupling between the magnetic moments localized at 
the edges can be controlled by means of either electron 
or hole doping of the nano ribbons ( Jung and MacDonaldl . 
120091; ISawada e* all 120091 ) . High doping levels eventually 
suppress magnetism since the flat band shifts away from 
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FIG. 10 (a) An atomically resolved STM image showing 
the presence of both armchair and zigzag graphene edges, 
(b) Typical STS curve measured at a zigzag edge reveals 
the peak related to ze ro-energy states. Reproduced from 
ijKobavashi et all I2005T ). Copyright 2005 by the American 
Physical Society. 



the Fermi level, thus eliminating the electronic instabil- 
i ty associated with the prese nce of low-energy electrons 
(|Jung and MacDonaldl 120091) . 

It is worth mentioning that at the time this review was 
written, no direct proof of edge magnetism in graphene 
has been reported. However, the presence of local- 
ized low-energy states at zigzag edges of graphene has 
been verified by means of scan ning tunneling microscopy 
IjKobavashi et all I2006L l20o1ih . 



B. Possible applications in spintronics 

It has been realized that the intriguing magnetic prop- 
erties of graphene nanostructures may find applications 
in spintronics. The pioneering idea was introduced by 
Son, Cohen and Louie, who have predicted that external 
electric fields induce hal f -metall icity in zigzag graphene 
nanoribbons ( Son et all , [2006b). The half-metallicity 
refers to the coexistence of a metallic state for electrons 
with one spin orientation and an insulating state for elec- 
trons with the opposite spin orientation. An electric 
field is applied across the nanoribbon as shown in Fig- 
urc fTTT a). At zero field the system is characterized by the 
energy gap A5 for the spin-polarized states localized at 
both edges (Fig. fTTTb')'). An applied electric field breaks 
the symmetry and closes the gap for one of the spin direc- 
tions selectively (Fig. [TTJc)). The critical field required 




e V>0 



FIG. 11 Scheme of electric-field-induced half-metallicity in 
zigzag graphene nanoribbons. (a) Electric field is applied 
across the nanoribbon, from left edge (spin-up, red arrows) 
to right edge (spin-down, blue arrows), (b) Schematic repre- 
sentation of the spin-resolved local density of states for the 
opposite edges at zero applied field, (c) Applied electric field 
closes the band gap at the Fermi level Ef for spin-down elec- 
trons selectively. 
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FIG. 12 Spin injection from (a) an ideal zigzag graphene 
nanoribbon and (b) a nano ribbon with a defect ive bottom 
edge into n-doped graphene (|Wimmer et al. 1. 120081 ). Nonequi- 
librium densities for spin-up and spin down-electrons are 
shown in blue and red, respectively. Copyright 2008 by the 
American Physical Society. 



for inducing the half-metallicity is 3.0/u> Volts, where w 
is the nanoribbon width in A. The direction of the ap- 
plied electric field defines the spin channel with metallic 
conductivity. If realized in practice, this simple device 
would offer efficient electrical control of spin transport - 
a highly desirable component for spintronics. 

In addition to the device described, several other ap- 
proaches for controlling the spin transport in graphene 
nanostructures have been proposed. One of them ex- 
ploits disorder for ach ieving the goal; an example from 
(jWimmer eiaJl l200Sh is shown in Figure [12] Electric 
current flowing along the edges of zigzag nanoribbon in- 
jects spin-polarized electrons into a graphene reservoir 
(Fig. HHa)). However, the net spin-polarization of the 
current is zero due to the antifcrromagnetic coupling be- 
tween the two equivalent edges. Then, an extended de- 
fect is introduced into one of the edges as shown in Fig- 
ure [T27 b) . The defect both quenches magnetic moments 
and scatters the carriers at the rough edge. However, 
conduction at the opposite edge remains unaffected, thus 
allowing for injecting a current with a net spin polariza- 
tion. Other proposals based on defect and impurity engi- 
neering of spi n transport in graphene nanoribbon s have 
been reported (ICantele et al\ l2009tlLakshmi et all 120091: 
iPark et all l2009t iRocha et all |2009|) ~ 

Magnetic graphene nanostructures were also proposed 
as components of magnctorcsistive junctions. Such de- 
vices are currently used as magnetic field sensors, e.g. 
in the read heads of hard disk drives. Typical magnc- 
torcsistive junctions involve ferromagnetic metal layers 
separated by a nonmagnetic spacer layers, e.g. Co layers 
separated by a non-magnetic Cu layer, or bec Fe layers 
separated by a few-nanometers thick layer of crystalline 
MgO. A crucial characteristic of such spintronic devices 
is their magnetoresistance ratio (MR) which shows the 
change in electric resistance as a function of the rela- 
tive orientat ion of the magnetiz ation of two ferromag- 
netic layers ( Heiliger et all [2006h . This quantity can be 
defined as 





MR: 



R AP - R p 
min(i? p ,i? AP ) 



x 100%, 



(9) 



FIG. 13 Scheme of magnetoresistive device based on zigzag 
graphene nanoribbon connecting two ferromagnetic contacts 
(jKim and Kiml . l2008h . The low-resistance configuration (a) 
and the high-resistance configuration (b) of the device corre- 
spond to parallel and antiparallel orientations of the magnetic 
moments of the two ferromagnetic leads, respectively. Arrows 
denote the magnetic moments of both graphene edges and fer- 
romagnetic leads. 



where R p and R AP are the resistances for parallel and an- 
tiparallel relative orientations of the magnetic moments 
of the layers. Magnetoresistive devices with high mag- 
nitudes of MR are demanded by the future nanoscale 
electronics. It has been predicted that a zigzag graphene 
nanoribbon placed between two ferromagnetic contacts 
constitutes a magnetoresistive ju nction with very hig h 
values of magnetoresistance ratio ( Kim and Kiml . ]2008[) . 
The low-resistance state of the proposed device corre- 
sponds to the parallel configuration in which the mag- 
netic moments at graphene edges are coupled ferromag- 
netically to each other. The ferromagnetic coupling is 
enforced by the strong interaction with the magnetic 
moments of ferromagnetic contacts (Fig. Eta)). In 
the antiparallel configuration (Fig. fT3Tb)). the magnetic 
graphene nanoribbon develops a domain-wall arrange- 
ment of edge spins with high resistance. It is worth men- 
tioning that spin-transport measurements in micrometer- 
scale lateral graphene devices contacted by ferromag- 
netic electrodes have been carried out experimentally 
dHill et all 120061 : 1 Jozsa et all 120081 : pTombros et q"H . l2007l 
l2008h . However, the magnetoresistance effect observed in 
these experiments is due to the long spin-diffusion lengths 
in graphene. 

An all-graphene device based on the armchair-zigzag- 
armchair nanoribbon junction has also been predicted 
to sh ow magnetoresistance effect (jMuiioz-Rojas et all 
2009). In zero applied magnetic field the magnetic cou- 
pling between the opposite zigzag edge segments is an- 
tifcrromagnetic (Fig. EJa)) and, hence, electric resis- 
tance is high due to the gapped electronic state. Suf- 
ficiently strong magnetic fields favor the parallel config- 
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(|Reina et aU l2009h . 



Considerable progress has also been achieved in con- 
trolled manufacturing of graphcne nanostructures. The 
lithographic patterning allows to produce graphen e 
nanoribbons as narrow as ^15 nm (lHan et all l2007h . 



Sub-10-nanometer graphene nanoribbons have been syn- 
thesized using a vari ety of chemical approaches starting 
from either graphite (iLi et all [2001 IWang et all l2008bh 
or ca rbon nanotubes (jjiao et all 120091 : iKosvnkin et all 
I2009D . Well-ordered edges along a single crystallographic 
direction of graphene have be en produced by means 



of ch emical vapor deposition ( Campos-Delgado et all 
120081 ). annealing by Joule heating ( Jia et aL . l2009l ). and 
anisot ropic etching of grap hene using metallic nanopar- 
ticles ( Campos et all I2009T) . 




FI G. 14 AU-graphene maenetoresistive device proposed 
in (|Munoz-Roias et all [2009). The device represents an 
armchair-zigzag-armchair nanoribbon junction. The high- 
resistance antiparallel configuration (a) switches into a low- 
resistance parallel configuration (b) in a strong enough ap- 
plied magnetic field. 



uration (Fig. Qjjb)) which shows a lower resistance due 
to the spin-polarized edge states crossing the Fermi level 
(Fig. HJc)). This ultrasmall device thus acts a magnetic 
field sensor capable to detect magnetic fields from few 
hundreds of Gauss to several Tesla at low temperatures. 

Epitaxial graphene and the isostructural hexagonal 
boron nitride (/i-BN) have been proposed as efficient ul- 
trathin non-magnetic spacer s for traditional mu l tilayer 
magnetoresistive junctions dKaroan et all I2007L 120081 : 
lYazvev and Pasquarellol l2009h . Unlike in the device 
shown in Fig. II 31 mono- or multilayers of graphene or 
/i-BN are sandwiched between two ferromagnetic lay- 
ers. The transport direction is orthogonal to the plane 
of spacer layers. The key to feasibility of such de- 
vices is the fact that the lattice constants of graphene 
and /i-BN match closely those of Co and Ni. More- 
over, in the case of multilayer graphene the momen- 
tum selection criteria allow efficient transport only for 
the minority-spin channel in the parallel configuration 
of the device. Very high magnetoresistance ratios have 
been pre dicted for mult i layer graphe ne used as a spacer 
material (jKarpan et all , l2007U2008r) . Meanwhile, high- 
quality epitaxial monolayers of both graphene and /i-BN 
on ferromagnetic transition metals have been grown ex- 
perimen tally using the chemic al vapor deposition tech- 
niques (iDedkov et ai , 2~008al lbl: iGruneis and Vvalikh 



20081: lOshima and Nagashimal Il997t iRader et ail 12009 : 



Varvkhalov and Raderl . l2009UVarvkhalov et all L'OO.sli 
Successful growth of lower-quality multilayer graphene 
on polycrystalline Ni substrates has also been reported 



C. Magnetic ordering at finite temperatures 

R has already been mentioned at the beginning 
that the Curie temperatures of ferromagnetic materi- 
als must be higher than the operation temperature of 
the device, which is typically supposed to be close 
to 300 K. When introducing the working principles 
of the proposed spintronic devices based on graphene 
edges, temperature-related limitations were not dis- 
cussed. However, magnetic order in low-dimensional sys- 
tems is particularly sensible to thermal fluctuations. In 
particular, the Mermin- Wagner theorem excludes long- 
range order in one-dimensional magnetic systems (such 
as th e magnetic graphene edges ) at any finite tempera- 
ture ( Mermin and Wagnerl . Il966l) . The range of magnetic 
order is limited by the temperature-dependent spin cor- 
relation lengths £ Q (a = x, y, z) which define the decay 
law of the spin correlation 



) = (SfOexpH/D, 



(10) 



where sf are the components of magnetic moment unit 
vector Si at site i. In principle, the spin correlation length 
imposes the limitations on the device dimensions. In 
order to establish this parameter one has to determine 
the energetics of the spin fluctuations contributing to the 
breakdown of the ordered ground-state configuration. 

The energetics of the transverse and longitu- 
dinal spin excitations (Figs. [T5Ta).(b)) have been 
explored using density - funct ional-thcory calculations 
(lYazvev and Katsnelsonl . 120081 ). The magnetic correla- 
tion parameters in the presence of spin-wave fluctua- 
tions, the dominant type of spin disorder in this case, 
were obtained with the help of one-dimensional Hciscn- 
berg model Hamiltonian 



i+li 



(11) 



where the Hciscnbcrg coupling a = 1nja\ = 105 meV 
corresponds to the spin-wave stiffness k = 320 meV A 2 
calculated from first principles. The estimated small 
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reaching this goal include chemical functionalization of 
the edges with heavy-element functional groups or cou- 
pling graphcnc to a substrate. 
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FIG. 15 Schematic representation of the transverse (a) and 
the longitudinal (b) low-energy spin excitation at zigzag 
graphene edges. The magnetic moments are shown by arrows. 
Their directions are represented by orientation and color of 
the arrows, and the magnitude is shown by their lengths and 
color intensity, (c) Correlation lengths of magnetization vec- 
tor components orthogonal (£ z ) and parallel (£ x , £ y ) to the 
graph ene plane as a function of te mperature. Reproduced 
from (|Yazvev and Katsnelsonl . 120081 ). Copyright 2008 by the 
American Physical Society. 



anisotropy parameter d/a w 10~ 4 originates from the 
weak spin-orbit interaction in carbon. This simple mode l 
Hamiltonian has known analytic solutions ( Jovcd . Il967h . 
FigurcfTSTc) shows the spin correlation lengths calculated 
for our particular case. Above the crossover temperature 
T x pa 10 K, weak magnetic anisotropy does not play any 
role and the spin correlation length £ oc T . However, 
below T x the spin correlation length grows exponentially 
with decreasing temperature. At T = 300 K the spin 
correlation length (pal nm. 

From a practical point of view, this means that the 
dimensions of spintronic devices based on the magnetic 
zigzag edges of graphene and operating at normal tem- 
perature conditions are limited to several nanometers. 
At present, such dimensions are very difficult to achieve, 
which can be regarded as a pessimistic conclusion. Nev- 
ertheless, one has to keep in mind that the spin stiffness 
predicted for the magnetic graphene edges is still higher 
than the typical values for traditional magnetic mate- 
rials. That is, graphene outperforms d-element based 
magnetic materials, and there is a room for improvement. 
Achieving control over the magnetic anisotropy d/a could 
possibly raise the crossover temperature T x above 300 K 
and thus significantly extend £. Possible approaches for 



VI. MAGNETISM IN GRAPHENE AND GRAPHITE 

A. Radiation damage and defects in carbon materials 

Experimental observations of ferromagnetic ordering 
in irradiated graphite have already been mentioned in 
the introductory part of this review. These results are 
particularly exciting because of the fact that the induced 
magnetic ordering is stable at room temperature and well 
above. Let us now try to understand the origin of mag- 
netism in irradiated graphite. The present section cov- 
ers the cases of both graphcnc and graphite which has a 
three-dimensional crystalline lattice composed of weakly 
coupled graphene layers. 

The basic picture of the radiation damage process 
in carbon materials is relatively simple. Irradiation of 
graphite with high-energy particles (e.g. protons) pro- 
duces several types of point defects. In carbon materials 
the defects are created as a result of so-called "knock-on 
collisions" (|Banhartl[l999H K rasheninnikov and Banhart] . 
2007). This process involves the direct transfer of ki- 
netic energy from the high-energy incident particles to 
the individual atoms in material's lattice. If the trans- 
ferred energy is larger than the displacement threshold 
Td, the recoil atom may leave its equilibrium position 
leading to the formation of a pair of point defects - 
a vacancy defect and an interstitial. The structure of 
the vacancy defect in graphene and graphite is shown 
in Figure [Tffiah In gr aphene the interstiti al defects 
have a bridge structure ( Lehtinen et all l2003h . while in 
graphite the stable configuration corresponds to a car- 
bon atom trapped between the adjac ent graphene lay- 
ers (|Li et all 120051 ; iTelling et all |2003() as shown in Fig- 
ures Urir b). The displacement threshold Td for carbon 
atoms in graphitic materials was found to be ~20 eV in a 
number of studies ( Crespi et all 11996k ISmith and Luzzl . 




vacancy 



interstitial 



hydrogen 
chemisorption 



FIG. 16 Atomic structures of primary point defects pro- 
duced upon irradiating graphite by high-energy protons: (a) 
single-atom vacancy, (b) interstitial bridging the neighboring 
graphene layers, and (c) hydrogen chemisorption. The dan- 
gling bond (DB) of the reconstructed vacancy defect and the 
chemisorbed hydrogen atom (H) are labeled. 
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120011: lYazvev et all l2007t IZobelli et all 120071) . Cre- 
ation of defects due to electron stopping, i.e. the pro- 
cess involving possible electronic excitations and ion- 
ization of individual atoms, is less important in car- 
bon materials since electronic excitat ions in metals 
are delocalized and quench instan tly ( Banharti Il999t 
IKrasheninnikov and Banharti 120071 ). 

After slowing down, reactive particles may also pro- 
duce chemisorption defects. In particular, protons are 
able to bind to individual carbon atoms in graphene lat- 
tice resulting in their rehybridization into the sp 3 -state 
(Fig. ITST c)). Such defects are referred to as hydrogen 
chemisorption defects. From the point of view of one- 
orbital models that we use in our review, both vacancy 
and hydrogen chemisorption defects are equivalent. In 
both cases a defect removes one p z -orbital from the ir- 
system of graphene. In the first case, the ^-orbital is 
eliminated together with the knocked-out carbon atom. 
The hydrogen chemisorption does not remove the carbon 
atom from the crystalline lattice, but once rehybridized 
the atom is unable to contribute its p 2 -orbital to the tt- 
electron system. These two types of defects are further 
referred to as -vacancies. 

The defects described above are the primary defects 
in the radiation-damage process. More complex defects 
can be produced at later stages of the process. For in- 
stance, single-atom vacancies and interstitials may aggre- 
gate producing extended defects. Complexes involving 
two or more different defects can also be formed upon 
irradiation. Examples are complexes of hydrogen with 



irradiation. Hixampies are complexes or n;v 
vacancies and interstitials dLeht inen et all 



... 12004), and 

intim ate Frcnkel pairs (|Ewels et all l2003t lYazvev et all 
I2007D . Radiation damage in graphi tic materials may also 



produce the Stone- Wa l es de fects ( Kaxiras and Pandevl 
119881 : IStone and Walesl . Il986h . 



This resu lt has been wid e ly con f irmed using both first- 
principles (iDuplock et all 12004 iLehtinen efall . 12004 
lYazvev and Helm L 2007T) and mean-field Hubbard mode l 
( Kumazaki and Hirashimal 120071 : |Palacios~e£ 

Mi l2008h 

calculations. Figure IT7T a) shows the spin-resolved 
density-of-states plots for hydrogen chemisorption and 
vacan cy defects obtained using first-principles calcula- 
tions (jYazvev and Helml . 120071) . In the first case, the 
sharp peak close to the Fermi level corresponds to the 
quasi-localized state induced by the chemisorbed hydro- 
gen atom. The peak is fully split by exchange and the sys- 
tem is characterized by a magnetic moment of 1/^b at any 
defect concentration. The distribution of spin density 
around the defective site clearly shows a vo x \/3R30 o 
superstructure (Fig. IT7T b)). The case of vacancy defect 
is somewhat more complicated. In addition to the quasi- 
localized state, there is also a localized non-bonding state 
due to the presence of a cr-symmetry dangling bond in 
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B. Defect-induced magnetism in graphene 

The single atom p z -vacancies described above have a 
particularly profound effect on the electronic structure of 
ideal graphene. Let us consider a periodically repeated 
supercell of graphene composed of 2N (Na = Nb = N) 
carbon atoms. Elimination of one atom from sublattice A 
introduces a zero-energy state in the complementary sub- 
lattice (a = N B ; thus r) = 2N B - ((N A - 1) + N B ) = 1). 
Such zero-energy states extending over large distances 
are called quasi-localized state s since they show a p ower- 
law decay (jHuang et all . 120091 : iPereira et all . 120061 ) . The 
quasi-localized states have been observed in a large num- 
ber of scanning tunneling microscopy (STM) studies of 
graphite as triangular y/2> x V3-R30 superstructures ex- 
tending over a few nanom e ters and localized around point 
defects dKellv and Halai , Il998t iMizes and Fosterl . Il989t 



iRufReux et at . 2000h . For the single-defect model we 



have adopted, Lieb's theorem predicts a magnetic mo- 
ment of | (Na — 1) + N B \ = 1/Ab per supercell, that is, the 
presence of a defect induces ferromagnetic ordering. 




FIG. 17 (a) Spin-resolved density of states (DOS) for the 
vacancy and the hydrogen chemisorption defect in graphene 
calculated from first principles. Red and blue curves cor- 
respond to the majority and minority spins, respectively. 
Dashed curve shows the reference density of states of the 
ideal graphene. Zero energy corresponds to the Fermi level. 
Exchange-split peaks which correspond to quasi-localized 
(QL) and dangling-bond (DB) states are labeled, (b) Iso- 
surface representation of the spin-density distribution at the 
hydrogen chemisorption defect obtained from first princi- 
ples. Red and blue surfaces correspond to the majority- and 
minority-spin densities, respectively. 
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this defect fFig.fTBta)). The dangling-bond state shows a 
very strong exchange splitting and contributes lfiB to the 
total magnetic moment of the defect (Fig. 117( a)). How- 
ever, the magnetic moment due to the quasi-localized 
state is partially suppressed in this case due to the self- 
doping effect related to the structura l reconstruction of 
the vacancy (|Yazvev and Helml . l2007t ). The overall mag- 
netic moment per vacancy defect varies from 1.12/zb to 
1.53/xb for defect concentrations ranging from 20% to 
0.5%. 

Magnetic moments due to dangling bonds can 
also be contributed by other types of defects, e.g. 
the bridge-configurati on interstitial defect in graphene 
( Lehtinen et all |2003|) . However, one has to keep in 
mind that magnetic ordering due to only localized mag- 
netic moments in graphene-based system is improbable 
at high temperatures. The Ruderman-Kittel-Kasuya- 
Yoshida interaction is expected to be weak in this 
case due to the semi-met a llic e l ectronic structure of 
graphene (iBrev and Fertid. 20061 : Dugaev et all 120061 : 
ISaremil . [2007 : IVozmediano et all 12005( 1 . On the other 
hand, magnetic ordering due to the quasi-localized states 
can be considered as itinerant magnetism without exclud- 
ing a possible contribution of dangling-bond magnetic 
moments to the net magnetic moment of a defective car- 
bon system. 

The system with one defect placed in a periodically 
repeated supercell is only a rough model of disordered 
graphene for two reasons. First, all defects arc located 
in the same sublattice of the graphene layer. Second, the 
defects form an ordered periodic superlattice. A more 
realistic description of disorder can be achieved by con- 
structing models with de fects randomly distributed in 
a large enough supercell ( Yazvevl . l2008bl l. Such mod- 
els allow defects to occupy both sublattices at arbitrary 
concentrations and eliminate any short-range order in 
the spatial arrangement of defects. Larger supercells are 
needed for building the disordered models which makes 
first-principles calculations impractical. However, such 
system can still be treated using the mean-field Hubbard 
model calculations. 

Figure I18f a) shows the distribution of spin density in 
a selected region of a large supercell randomly popu- 
lated by p z -vacancies. Defects in different sublattices 
are shown as black triangles of different orientations. 
The resulting picture can be explained if one considers 
the following two arguments. First, from Lieb's theo- 
rem the total magnetic per supercell is M = \Na — 
n b\ = \N$ - N$\, where and N% are the num- 
bers of defects created in sublattices A and B, respec- 
tively. This means that electron spins populating the 
quasi-localized states in the same sublattice are oriented 
parallel to each other while the antiparallel arrangement 
is realized when electron spins populate different sub- 
lattices. This conclusion is qualitatively the same as 
the one we obtained when considering finite graphene 
nanofragments. Second, quasi-localized states populat- 
ing complementary sublattices interact with each other 
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FIG. 18 (a) Distribution of local magnetic moments Mi in 
a small region a system with random distribution of p z - 
vacancies obtained thought the mean-field Hubbard model 
calculations. Positions of defects belonging to sublattices A 
and B are marked by ▲ and ▼, respectively, (b) Average 
magnetic moments for the atoms in sublattice A (V,T) and 
sublattice B (A, A) as a function of defect concentration x. 
The defects are either distributed equally between the two 
sublattices (solid curves/filled symbols) or belong to sublat- 
tice B only (dashed curves/empty symbols). Net magnetic 
moments per carbon atom (dashed line) is shown for the case 
of defects distributed over sublattice B only. 



( Kumazaki and Hirashimal l2007t iPalacios et al. . 2008). 
The interaction lifts the degeneracy leading to weakly 
bonding and anti-bonding states. This provides another 
mechanism for escaping the instability associated with 
the presence of low-energy electronic states. The in- 
teraction strength between two defects increases with 
decreasing distance between them. For very short dis- 
tances, the gain in exchange energy does not compen- 
sate for the kinetic energy penalty due to the splitting. 
This leads to quenching of the de f ect-induced ma gnetic 
moments (|Boukhvalov et all [20081 : 1 Yazvevl . l2008bl ). The 
competition between these effects is demonstrated in Fig- 
ure QUa). 

More quantitative results are presented in Figure HHTb) 
which shows the mean magnetic moment (M A ) and 
(M B ) per carbon atom in sublattice A and sublattice B 
as a function of defect concentration x ( Yazvevl . 12008 hi ) . 
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The resulting values have been averaged over many ran- 
dom placements of defects in the simulation supercell. 
The plot refers to the situation of defects equally dis- 
tributed over the two sublattices (iVj[ = Ng) and to 
the situation when defects belong to sublattice B only 
(N% = 0). In the first case, the magnetic moments in the 
two sublattices compensate each other. The overall mag- 
netic ordering is of antiferromagnetic character. When 
defects populate only one sublattice the system exhibits 
ferromagnetic ordering. The net magnetic moment per 
carbon atom (M) = {{M A ) + (M B ))/2 = x/2 scales 
linearly with the defect concentration. Both numerical 
results are in full agreement with Lieb's theorem. 



C. Magnetism in graphite and multilayer graphene 

Experimental observations of ferromagnetism in irra- 
diated graphite point to the conclusion that sublattices 
of individual graphene layers in bulk graphite are pop- 
ulated by defects differently: that is, there must be a 
mechanism which makes the sublattices of graphene in- 
equivalent. Such an intrinsic discriminating mechanism 
was ascribed t o the stacking order of graphite layers in 
bulk graphite (|Yazvevl . l2008b|) . The lowest-energy ABA 
stacking order of individual graphene sheets in graphite 
breaks the equivalence of the two sublattices as shown in 
Figure H9T a). In fact, only local ABA stacking order is 
required to discriminate between the two sublattices of 
the middle sheet. The mechanism can be demonstrated 
for the case of hydrogen chemisorption defects. First- 
principles calculations show that the configuration which 
involves hydrogen chemisorbed on sublattice B is 0.16 eV 
lower in energy than on sublattice A (see Figs.fTWb) and 
[IW c)). This energy difference is sufficient to trigger a 
considerable difference in equilibrium populations of the 
two sublattices. The energy barrier for the hopping of 
chemisorbed hydrogen atoms is relatively small (~1 eV) 
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FIG. 19 (a) Inequivalent carbon atoms (A and B) in 
ABA graphite, (b) Possible pathways for the diffusion of 
chemisorbed hydrogen in graphite, (c) Schematic representa- 
tion of the potential energy surface for the in-plane diffusion of 
hydrogen in graphite showing the relative energies fo r the lo- 
cal mi nima and transition states. Reproduced from (jYazvevl 
l2008bh . Copyright 2008 by the American Physical Society. 



to allow for thermally activated diffusion at mild temper- 
atures. 

Similar discriminating mechanisms may also exist for 
the other types of defects created by irradiation, e.g. for 
vacancies. Cross-sections for momentum transfer dur- 
ing knock-on collisions with high-energy incident parti- 
cles are likely to be very similar for both A and B car- 
bon atoms in graphite. However, the stacking order may 
have a strong influence on the recombination of intersti- 
tial and vacancy defects close to equilibrium conditions. 
It was also shown that instantaneous recombination of 
low-energy recoil atoms in graphit e is significant l y mor e 
probable for atoms in position A ( Yazvev et all l2007t ). 
That is, more vacancies in sublattice B are left assum- 
ing an equal number of knock-on collisions involving the 
atoms of both types. These results allow us to conclude 
that the most probable physical picture of magnetic order 
in irradiated graphite is ferrimagnetism. The magnetic 
moment induced by defects in sublattice A is larger than 
the one induced in sublattice B. 

It is worth mentioning other possible scenarios for the 
onset of magnetism in graphene-related materials and 
nanostructures. It was shown that local negative Gauss 
curvatures in g raphene layers als o lead to localized mag- 
netic moments ( Park e^all . 12003). Coupling between the 
two graphene layers in a biased graphene bilayer intro- 
duces a pair of sharp van Hove singularities close to the 
Fermi level. It was shown theoretically that in a range 
of conditions biased bilayer graphene satisfies the Stoncr 
criterion lead i ng to a low-density ferromagnetic phase 
(|Castro et ^.1 . 120081 ). 

The question of magnetic ordering in defective 
graphene and graphite at finite temperatures remains 
largely unaddrcsscd. Similarly to the one-dimensional 
system discussed above, an isotropic two-dimensional 
system cannot develop lo ng-range magnetic orderin g 
at any finite temperature ( Mermin and Wagnerl . Il966l) . 
However, the introduction of a small magnetic anisotropy 
d/a ~ 1 0~ 3 already leads to ve r y high transition temper- 
atures ( Barzola-Quiauia et all 120071. Weak magnetic 
coupling between the individual layers in graphite also 
produces a pr onounced effect on the magnetic transition 
temperature ( Pisani et all 120081) . 



VII. CONCLUSIONS AND PERSPECTIVES 

The review illustrated a rich variety of magnetism sce- 
narios in graphene nanostructures and explained them in 
terms of simple physical models. Beyond these theoret- 
ical considerations the field of carbon-based magnetism 
faces a number of challenges. The most important prob- 
lems are related to the experimental side of the field. In 
particular, the physics of magnetic graphene edges has 
already attracted a large number of computational and 
theoretical researchers. However, no direct experimen- 
tal evidence has been reported at the time this review 
was written. Further progress in this field will also re- 
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quire novel manufacturing techniques which would allow 
control of the edge configuration with truly atomic preci- 
sion. The area of defect-induced magnetism in graphite 
demands detailed studies of defects produced upon irra- 
diation and their role in the onset of ferromagnetic or- 
dering. The limits of saturation magnetization and Curie 
temperature in irradiated graphite have still to be estab- 
lished. On the theory side of this field, an understanding 
of magnetic phase transitions in graphene materials and 
nanostructures has to be developed. Other important 
directions of theoretical research include spin transport 
and magnetic anisotropy of carbon-based systems. 
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